Signature of chirality in scanning-probe imaging of charge flow in graphene 
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We theoretically propose to directly observe the chiral nature of charge carriers in graphene mono- 
and bilayers within a controlled scattering experiment. The charge located on a capacitively coupled 
scanning probe microscope (SPM) tip acts as a scattering center with controllable position on the 
graphene sheet. Unambiguous features from the chirality of the particles in single and bilayer 
graphene arise in the ballistic transport in the presence of such a scattering center. To model 
the scattering from the smooth potential created by the SPM tip, we derive the space-dependent 
electron Green function in graphene and solve the scattering problem within the first-order Born 
approximation. We calculate the current through a device with a SPM tip between two constrictions 
(quantum point contacts) as a function of the tip position. 
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The isolation of few and single layer graphene,i^i^ the 
two-dimensional carbon allotrope, triggered tremendous 
research activities (for a review, see Ref. |3) . Graphene 
is technologically of high interest f^i^ii as it is chemically 
stable, but is also very appealing for fundamental sci- 
entific research. The electronic eigenstates in graphene 
obey a linear gapless dispersion relation around the Fermi 
energy, and bear a pseudospin, which is always aligned 
parallel in the direction of momentum. This properties 
imitate the behavior of chiral masslcs Dirac particles^^i^ 
discussed in relativistic quantum mechanics. Therefore 
many concepts of solid-state physics are now being recon- 
sidered for pseudo-rclativistic carriers while at the same 
time, effects known from relativistic quantum mechanics 
can be found in solid-state physics. Examples for this 
are the unusual energies of the Landau levels and the 
Klein paradox^ Moreover, long spin relaxation lengths^ 
make graphene an interesting system for spintronicsi^ 
and spin-based quantum information processingJ^ 

The high carrier mobility of graphene has lead to 
an active discussion of impurity scattering. Cur- 
rently, it is assumed that scattering from Coulomb 
potentialsi^ii^iiiii^ii^iiiii^ii^i^ai^i limits the conduction 
electron mobility, while short-ranged defects are less 
relevant 1^ 

In the following, we discuss the possibility for a con- 
trolled experiment to test whether the charge carriers 
in graphene behave like chiral particles in a scattering 
event. We propose to use the method of mapping electron 
flow by scanning-probe microscopy (SPM) as developed 
and applied to two-dimensional electron gases (2DEGs) 
in semiconductors by the Westervelt group.— Topinka 
et alM^ demonstrated that coherent electron flow in a 
2DEG formed in a GaAs/GaAlAs heterostructure can be 
imaged directly by placing a charged SPM tip on top of 
the sample. The tip, being capacitively coupled to the 
sample, repels the conduction electrons beneath, forming 
a circular scatterer with a precisely controllable position. 
By scanning the tip over the conductor, the conductance 
of the sample is modified depending on the current den- 
sity beneath. By setting the SPM tip directly behind 
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FIG. 1: The proposed setup. An applied bias voltage V be- 
tween the source S and the middle region M injects a current 
into the graphene sheet. The injected electrons scatter from 
an artificial scatterer created by a SPM tip above the surface. 
The flux of electrons that are coherently scattered into the 
drain D can be detected in the drain current I. 



a constriction (quantum point contact, QPC),— the 
dominant mechanism of the conductance change is di- 
rect backscattering through the QPC next to the source, 
see Fig. [1] In single layer graphene, this backscattering 
is forbidden;^ and therefor one can not expect any re- 
sistance change between the respective contacts S and 
M in such a setting. To actually use forward scattering 
from the tip, one might consider putting the tip in front 
of the constriction. Such a setup, however, suffers from 
the uncontrolled direction of propagation of the incom- 
ing particles. To control the direction of propagation of 
in and outgoing particles, i.e., the scattering angle, we 
propose to use two constrictions (QPCs), as in Ref. [2^. 
By applying a voltage V between regions S and M, a cur- 
rent is injected into the middle region (M). The coherent 
scattering from one QPC to the other gives rise to a mea- 
surable current I in the drain D, which can be precisely 
measured as a voltage between M and The middle 
region M is supposed to be large, acting as a reservoir, 
absorbing all electrons that are not scattered into the 
drain D. Such an experiment could directly probe the 
differential scattering cross section for pseudorelativistic 
chiral particles and thus it demonstrates the chiral nature 
of the particles. 
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FIG. 2: Formulation of the scattering problem: An incoming 
plane wave propagating along the x-direction hits the scatter- 
ing potential. A detector measures the outgoing flux of the 
scattered wave as function of the deflection angle 0. 



For the microscopic description of scattering, and in 
the closely related Coulomb impurity problem,— the 
method of partial-wave expansion was adapted.— i^i^i^ 
This method is primarily suitable for strong short-ranged 
potentials such as impurities. To describe the scattering 
from the weak potential created by a SPM tip. wc apply 
the method of first-order Born approximation by deriv- 
ing the Green functions for single and bilayer graphene, 
in real space representation!^ 

For single layer graphene, the free Dirac Hamiltonian 
for the envelope wave function at the if-point is given by 
H ~ —ihvpic^xdx + cTydy) and in cylindrical coordinates 
r = (r, (j)) by 
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The scattering solution of the Dirac equation [H + 
U{r)]ip = Ekip, with Ek = hvpk, can be constructed 
using the Green function 
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Ho(kp) —iH^i{kp)e 



, (2) 



where p = r — r' = {p,0). The Green function is a 
solution of {H — E)G = —5{p), satisfying the outgoing 
radiation condition; thus, the use of the n-th order Han- 
kel function Hn{z) = The functional form of the 

Green function is consistent with Ref.[3l|. Since the Dirac 
Hamiltonian is formally closely related to the Rashba 
spin-orbit interaction, the derivation of the Green func- 
tion in Eq. Q is analogous to Refs. [s^lssl . 

The idealized scattering experiment is depicted in 
Fig. [21 A chiral plane wave il^vtif) = e*'^^[l,l]"^ prop- 
agating along the (so chosen) x-axis hits the scatter- 
ing potential. After the interaction with the potential, 
a detector measures the flux of the scattered wave as 
a function of the deflection angle <(). Within first-order 
Born approximation, the total wave function 'ip{r) of an 
electron which scatters at a potential U{r) is given by 
il}(r) = ■(/)o(r) + j <Pr'G{r - r')U{r')ipo{r'). Far away 
from the scattering center the wave function has the 
asymptotic form ip{r) = ipo{r) + /((/>) e''^'' [1, e^'^]'^ / ^/r. 



The scattered wave describes a chiral circular outgoing 
wave. The scattering amplitude is given by 
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with the Fourier transform U{q) = J(Pr' C/(r')e"^'' . 
The orientation of the vector q — k{ek — e^) of the mo- 
mentum transfer during the scattering process is deter- 
mined by the unit vectors and ej, in r and k direction, 
while its magnitude is given by g = 2fc sin((/)/2). The 
functional form of the scattering amplitude is in close 
analogy to nonrelativistic two-dimensional scattering.— 
Two main differences appear. First, the prefactor 
g-i0/2 t leads to the absence of any backscattering 
at potentials, irrespective of the potential shapci^ Sec- 
ond, the forward scattering is enhanced by a factor of 2 
compared to classical electrons. In gapped 2DEGs, scat- 
tering at potentials smaller than the Fermi wave length 
is dominated by s-wave scattering. For chiral particles in 
graphene, the free angular momentum eigenstates bear 
half-integer angular momentum related to the appear- 
ance of a Berry phase of ttJ^ Therefore, the scattering 
amplitude is dominated by states with orbital angular 
momenta of -1-1/2 and —1/2. In backward direction, the 
interference of the two states is destructive, leading to a 
suppression of scattering. In forward direction, the inter- 
ference is constructive, enhancing the scattering ampli- 
tude by a factor of 2. For scatterers larger than the Fermi 
wavelength (beyond the range of validity of the Born ap- 
proximation) , one can expect that the enhancement of 
forward scattering becomes non-universal, depending on 
the details of the potential. 

By comparison with the exactly solvable problem of a 
stepwise constant potential, the range of validity for the 
Born approximation can be estimated as (fci?)^ < Ek/Uo, 
where R and Uq denote the characteristic potential size 
and strength. As an interesting side note, for the physi- 
cally relevant Coulomb potential, the first-order Born ap- 
proximation generates accidentally good results, in three 
as well as in two dimensionsi^ 

When we restrict the calculation to circularly symmet- 
ric potentials, the scattering amplitude Eq. ([3]) can be 
further simplified using 



C/(q) = 27r dr' r' Joiqr')U{r') 



(4) 



The integral in Eq. ^ can be solved analytically for a 
variety of different scattering potentials, including an un- 
screened and the exponentially screened Coulomb po- 
tential, the stepwise constant potential, potentials of 
Lorcntzian and Gaussian forms, and the potential of a 
point charge above the graphene sheet 

In graphene bilayers, arranged in an A-B (Bernal) 
stacking, the charge carriers behave like massive gap- 
less Dirac fermionsi^ To compute the Green func- 
tion, one can directly make an ansatz motivated by 



the observation that the operators a 
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:9<i) and 



-ie^'^idr + '-d^) in Eq. © act 



as a type of ladder operators on the Hankel func- 
tions, i.e., atir„(/cp)e™'' = +iA:i/„+i(fcp)e*("+i)« and 
aff„(fcp)e"^ = -«fci?„_i(fcp)e'("-i)^^ Note the rela- 
tive minus sign for and a. With this representation, 
the Hamiltonian for the bilayer as derived in Ref. [s^ is 
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and the basic functional structure of the Green function 
for particle energy Ek = ft^fc^/2m directly follows as 



H2{kp)e^'' Hoikp) 



(6) 



The overall prcfactor can be deduced usingiiS that 
V4(fc2 + V2)iJo(fcp) = ~Hp) and the fact that aa'' = 
a^a = — V^. Using the Green function Eq. ([6]), we now 
construct the scattering wave function. Within the first- 
order Born approximation, the wave function takes the 

asymptotic form tp{r) ~ Vo(r)+/('?!>)e*'""'" [l, — e^*"^] /\/r. 
The scattering amplitude /(0) is then given by Eq. ([3]) 
with the substitution 2(j). 

For a circularly symmetric potential, the scatter- 
ing amplitude can again be simplified further by using 
Eq. (|4]) . The cases of non-relativistic particles in a 2DEG 
(j = 0), and of single layer [j = 1) and bilayer (j = 2) 
graphene are distinguished by their respective factors 
g-u0/2 (.Qg('j^^2) which are due to the Berry phase jtt 
acquired during the adiabatic propagation along a closed 
orbit. For a Gaussian potential U{r) = Uoe~^ 1"^^ , the 
resulting cross section da/dcj) = becomes 



da 



(X 



g-[2i?fesin(,/./2)]2 ^^^g2 n 
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We plot this result in Fig. [3l While the Berry phase pro- 
hibits backscattering in single layer graphene, backscat- 
tering is allowed in bilayer graphene, while scattering 
by an angle of ±7r/2 is forbidden in bilayers. These 
results are equivalent to the calculation of the scatter- 
ing cross section via the k-dcpcndcnt Green function, sec 
Ref. mili. 

To describe the SPM experiment, as shown in Fig. [2 
we model the potential of the charged tip as a Gaussian. 
Moreover, we consider the QPC at the source S as point 
source of chiral electrons (as determined by the Green 
function) at position {x^y) = (0,0). In the proximity 
of the scatterer, we approximate the incoming spherical 
wave ■(/; (centered around the QPC) as a plane wave and 
derive the scattered outgoing spherical wave (centered 
around the SPM-tip). The drain current / is calculated 
as the current at the location of the drain QPC at (x, y) = 
[d, 0) , where d is the distance between the source and 
drain QPCs. The normal component of the current is 
given by Jx — vpip'^axil^ for single layer, and by = 
~—lm\ili'' [a^dx + aydy)^] for bilayer graphene. Here ij} 
labels the sum of the incoming and outgoing spherical 
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FIG. 3: Differential cross section normalized by 
27r(/c_R)*(f/o/i^)^ in units of 1/fc for a circular scatterer 
of Gaussian shape. For single layer graphene, the Berry 
phase of tt prohibits backscattering while for bilayer graphene, 
it the Berry phase of 2n prohibits rectangular deflection. 
The radius of the scatterer is chosen kR = 0.5 for this plot. 



waves. In Fig. 51 we plot the relative change of the drain 
current AI / Iq due to the presence of the tip as a function 
of the tip position {x,y). Here, AI = I — Iq where Iq is 
the current in the absence of the SPM tip. 

In the presence of the scatterer, two ballistic trajecto- 
ries lead from source to drain: cither the electrons travel 
directly form the source (S) to the drain (D) or they scat- 
ter at the SPM tip potential and from there into the drain 
(D). The spatial pattern due to the interference between 
these two trajectories reveals the Fermi electron wave- 
length Xf ~ 2'K/k = hvp/Ep, the degree of coherence, 
as well as the scattering phase4£ 

The required scattering angle to pass from source to 
drain is a function of the tip position relative to source 
and drain. Therefore this experiment realizes to some 
extent an angle-resolved measurement of da/dcj). We as- 
sume the middle region (M) of the graphene sample to be 
large (but not larger than the coherence length Z^), i.e., 
R d <^ l^, so that scattering events with other angles 
will not significantly contribute to the drain current. For 
single layer graphene, one can expect a rapid loss of sig- 
nal if (f) > 7r/2, as scattering by larger angles is strongly 
suppressed. For bilayer graphene, one can expect for- 
ward as well as backscattering. However, a scattering 
angle of (/) = 7r/2 is forbidden. Therefore, a circular line 
appears (according to Thales' theorem), indicating a to- 
tal absence of scattering. In contrast, for a conventional 
semiconductor-based 2DEG, the intensity distribution is 
much more homogeneous. 

As shown in recent experiments;^ the minimal tip- 
induced potential width is about 300 nm at a potential 
height of 1 mcV. At a Fermi energy of 10 mcV in the 
graphene sheet, the range of validity of the first-order 
Born approximation (kR)^ < Ek/UQ is violated by about 
one order of mangitude. This violation leads to devia- 
tions mainly for forward scattering. Since in the pro- 
posed experiment the large angle scattering is of primary 
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FIG. 4: Relative change AI/Iq of the drain (D) current as 
a function of the tip position {x,y) in units of Af- The ex- 
istence of two possible trajectories, one directly from source 
(S) to drain, the other via the scatterer at {x, y), generates an 
interference pattern. Since scattering under an angle of 7r/2 is 
forbidden for a bilayer, a circle without signal (dashed white 
line) appears. The scattering potential radius and strength 
were chosen to be Rk = 1, and Uo/E = 0.3 respectively. 

interest, first-order Born approximation is still expected 
to deliver qualitatively correct results. 

Furthermore, sample roughness and disorder will also 
significantly modify this idealized experimental result, as 
already observed in 2DEGs.— i^i^^ Even so, if one as- 
sumes that small-angle disorder scattering is dominant, 
then disorder will be of most importance, when the tip 
is directly in between the source and the drain. With 



increasing scattering angle, the experimental result will 
be more and more robust against weak disorder. 

Our calculation was done for one of the two degener- 
ate valleys (the Dirac point at momentum K). For the 
other valley {K'), only the sign of (f> must be reversed in 
Eq. Therefore, the results for the scattering cross sec- 
tion Eq. ([7]) and the current (Fig. |4]) remain unchanged 
for K', as they are even functions of 0. Therefore, we 
expect these results to persist for arbitrary incoherent 




FIG. 5: Alternative setup where the QPCs have been replaced 
by two additional SPM tips playing the role of the source and 
drain contacts. 



mixtures of K and K' without any loss of interference 
visibility. 

To test the ballistic current through a graphene sheet, 
it would also be possible to use a multi-tip setup, see 
Fig. [5l as developed recently4^ Thereby, the roles of the 
source and drain contacts are played by two additional 
SPM tips. The SPM tip which is capacitively coupled 
and creates the scattering potential can also be moved 
behind the other tips and thus allows for a mapping of 
the complete angular dependence of the scattering am- 
plitude. The small cross sections of the contacts on the 
graphene ensure the angular resolution of the scattering 
experiment. This alternative setup, however, bears the 
experimental disadvantage that the resistance change as 
function of the scatterer position can not be measured in 
a non-local geometry as in the case of the QPC-setup. 

In conclusion, we propose to test the chirality of elec- 
trons in graphene mono- and bilayers in potential scat- 
tering, probed by a SPM tip in a transport setting with 
two QPCs. We describe the scattering within first-order 
Born approximation, which requires the derivation of the 
electron Green function in real space. 
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